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Abstract 



We calculate g A and g A in the Chiral Quark-Soliton Model (or equiva- 
lently in semibosonized Nambu-Jona-Lasinio model) in the limit of small 
and large soliton size. In small soliton limit we recover the results of the 

Constituent Quark Model. The agreement between the two models is 

(3) 

achieved due to the recently calculated 1 /N c contributions to g)^ and to 

. In the case of g^' these terms arise from the time-ordering of the 
collective operators. For large solitons simple scaling arguments show 
that the 1/N C contributions vanish with the inverse power of the soliton 
size. 
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The Constituent Quark Model (CQM), although often considered as naive, has been 
often used to calculate a priori dynamical hadronic matrix elements in terms of simple 
group-theoretical factors. The most prominent prediction of the CQM, namely the 
value of the singlet axial constant g^ — 1, caused the so-called "spin crisis". The CQM 
prediction [|IJ for = (N c + 2)/3, which overshoots the experimental value of 1.26, 
triggered the interest in the possible renormalization of a quark axial coupling away 
from 1 []2]|| . On the other extreme Skyrme-like models, in which baryons are considered 
as nontrivial configurations of the classical pion-sigma chiral fields, gave numbers in 
flat contradiction with the CQM (and experiment): g^ < 1 1 and g^ = |J. These 
models are subjected to a semiclassical quantization of the rotational zero-modes - only 
then the baryon spectrum is reproduced Q. The class of interpolating models, like the 
quark Nambu-Jona-Lasinio (NJL) model (or equivalently Chiral Quark-Soliton Model) 
PHlT|, in which both quarks and chiral fields are present, seemingly did not improve 
upon the <?a's, being on the low side in rough agreement with the Skyrme model. 
However not only the values of the #a's have been wrong, but also the N c counting. In 
the NJL model and in the Skyrme model as well, g^ scales like N c , however the 0(1) 
terms, which in the CQM are of utmost importance, seemed to vanish identically. The 
same concerns the 0(1) contributions to g^ which in the CQM sacles like 1. 

Recently 0(1) corrections to the ^a's and to the magnetic moments have been 
calculated in the NJL model [|i^-|i"5[ and also in the chiral bag model [Tj| . Not only have 
they significantly improved the numerics of the NJL model, but also the N c scaling, 
as predicted by the CQM, has been recovered. On the contrary in the Skyrme-like 
models based on a local effective meson lagrangian these corrections vanish identically 
Ref. 

The origin of these - for a long time overlooked - corrections (as has been shown in 



detail in Refs. [13—15]) is twofold : 1) the NJL model is formulated in terms of a non- 
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local effective action, namely in terms of the fermion determinant, and the operators 
appearing in the path integral for matrix elements have to be time- ordered, 2) the 
seemingly commuting c-number quantities like angular velocities and the Euler angles 
in fact do not commute, since upon the semiclassical quantization they are promoted 
to quantum operators. Apart from these time-ordered corrections there are also terms 
which come from the anomalous (or imaginary) part of the effective action (as in the 
case of ) . In this order of the 1/N C expansion their emergence does not depend on 
the time order of the collective operators. 

In this short note we show, that in the small soliton limit of the NJL model the 
results of the CQM for the #a's are recovered indepently on the model parameters. 
Indeed, when the soliton ceases to exist the quarks, which are normally bound in a 
potential created by the chiral fields, reach the plane wave orbits of the free Dirac 
particles |§ . Moreover, the contribution of the Dirac sea (or equivalently of the chiral 
fields) vanishes in this limit. Strictly speaking the Constituent Quark Limit can be 
accomplished by tuning the constituent quark mass m — > 0, in such a way that for each 
m the selfconsistent soliton solution is obtained. Since for some critical m cr the bound 
state solution ceases to exist even if the soliton has finite size, we choose first to push 
ro — > artificially for some unspecified m > m cr , for which the selfconsistent soliton 
solution does exist. It is then enough to observe, that the result does not depend on 
m. 

Furthermore, on the basis of a simple scaling argument we show, that in the large 
soliton limit the new corrections die out. Therefore in this limit the predictions of local 
mesonic effective theories are parametrically reproduced. The picture which emerges is 
therefore very appealing: the small and large soliton limits of the NJL model coincide 
with the previously obtained results within the CQM and the Skyrme model. Were the 
new non-local 0(1) corrections equal to zero, the nice correspondence with the CQM 
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would be completely lost. 

Let us first briefly remind the philosophy and main ingredients of the NJL model. 



Our starting point is a semibosonized Euclidean action of the NJL model [13J15 



S eS [U) = -Splog{-z $ + mlP B }, (1) 

where 7 matrices are antihermitian, the SU(2) matrix U 15 = exp(i7 5 T7r /F n ) describes 
chiral fields and m is the constituent quark mass. Eq. (|1|) has to be regularized. For the 
hedgehog Ansatz Uq 5 = cos P(r) + inTjs sin P(r) with the profile function satisfying 
the boundary conditions .P(O) = it and P(oo) = 0, the model has solitonic solutions. 
The existence of the soliton is achieved due to the interplay between the valence level 
and the energy of the continuum states. All quantities have therefore typically two 
contributions which we subsequently call valence and sea respectively. 

If the soliton size r$ (defined through P = P(r/ro)) decreases, the valence level 
joins the upper continuum and the sea contribution vanishes. For large solitons the 
valence level sinks into the Dirac sea and the explicit contribution of the valence level 
disappears. In this way the model interpolates between the CQM and the Skyrme-like 
models []. 

The energy of the soliton is given as the regularized sum of the eigenenergies of the 
Dirac hamiltonian: 

H(U ) = -M-ijkd k + m(J?). (2) 

H commutes neither with the isospin (T), the spin (S), nor with the angular momentum 
(L), but only with the so called grand-spin G = T + J = T + L + S. 



1 One should stress, however, that the four-derivative term from the NJL model does not 
coincide with the one of the Skyrme model. It is cutoff sensitive and it can stabilize the 



soliton only for very large constituent masses [10]. 
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To appreciate the way the constituent quark limit of the NJL model is actually 
reached, requires a careful analysis of the wave functions of the hamiltonian @. To 



this end one defines spin-isospin spherical harmonics [18 



~(±,±) 



G,M;j = G±±l = G±l), 
G,M-J = G±~,l = G >, 



(3) 



where G, M, j and I denote the eigenvalues of the grand-spin and its third component 
and of the total momentum and the angular momentum respectively. The soliton Dirac 
equation splits naturally into two parity V = II(— ) G sectors. Natural parity II = (+) 
spinors are given by: 

?(-.+) \ 



iAq{t) ^g'm + iC G ~(r) eg ., 



B+(r) - (+ ' +) 



-'G.M 



Dq{t) 



(4) 



whereas for II 



they take the following form: 



^~gA x ) 



?(-,+) 



(5) 



Then, for II = (+), the Dirac equation is reduced to a system of 4 first order differential 
equations for the radial functions A . . . D: 
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(6) 



where b G = 2^G{G+ 1). Grand spin subscripts G,M and parity superscript "+" 
have been suppressed. For II = (— ) one should change the sign of m in (^J). For 
any soliton shape one can choose two different boundary conditions in the origin: 1) 



A 



B 



M+1 and C = D = 0ot2)A = B = and C 



,K 



D 



Therefore the solutions are labelled by a multiindex n = (n, G, M, II, i), where n stands 
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also for a radial quantum number and % corresponds to the initial conditions discussed 
above. 

In order to quantize the model one performs collective quantization in terms of the 
rotation matrix A(t): U^(x) = A(t)U? (x) (t) , where A^ dA/dt = i/2 ELi^ a 
[ffl . As a result one gets a collective hamiltonian: 

n coU = M sol + — . (7) 

Here M so \ is the soliton mass and / the soliton moment of inertia. Spin operators J 
are related to the angular velocities by the quantization prescription: IQ — > J. The 
eigenstates of (0) have been successfully interpreted as baryons. Their wave functions 
are given m terms of D^ 2 \A) Wi gner matrices and the moment of inertia I is given 
as the sum over the eigenstates of the intrinsic Dirac hamiltonian (0) (in contrast to 
the collective hamiltonian (0)). Only the valence part of / is here of interest: 

N c „ (val | f | n)(n | f | val) 

Let us note in passing that since / ~ iV c each power of Q counts as one power of 1/N C . 

Using the same path-integral formalism, which has been employed to derive 7i co ii 
one can derive a compact expression for any quark operator of interest. The space-like 
collective axial current operator can be expressed by the following functional trace in 
Euclidean space fT3|j : 

with = 747J75T& and D a b = 1/2 Tr ( A^ r a Arb). For a singlet axial current Tb should 
be replaced by a unit matrix and one should also replace D ab — > l ab in eq.(|9|). It is 
implicitly understood that the sea part of eq.(|^) is regularized and that the vacuum 
contribution, corresponding to Uq 5 — > 1 is subtracted. 

Let us expand eq.(|5]) in Q, remembering that each power of Q counts as 1/N C . In 
the zeroth order we get: 



= i tr {r5 D ab ^-j^ = D aj a. (10) 

Note that at this level the singlet axial current is equal to in agreement with the 
Skyrme model result ||. 

The next term in Q (or in 1/N C ) is of great importance. Expanding in Q we get: 

There are two subtleties connected with eq.(|ll|): first it should be remembered that 
because of the collective quantization Q c is no longer a c-number but rather an oper- 
ator J c /I which does not commute with D ab . Second, the trace in eq.(|TT|) should be 
understood as time-ordered [P"3}-|H> . 

To this end let us consider the Euclidean propagator: 

< x I a 1 ii \ y> = - ty) E ®n{x)& n {y) exp(-E n (t x - t y )) 
vt + tt En>0 

- 0(t y - t x ) $n(£)<(0) exp{-E n (t x - t y )). (12) 

E n <0 



If two such propagators are multiplied and time ordering in eq.(|TTl) is assumed, then the 
correct expression for the baryon number one solutions reads (back in the Minkowskian 
space) : 



( A J) (1) = ^ E Y ^{<k\r b J \n><n\r c \k>D ab J c 

-En>M 

+ < n | V) | k >< k | r c | n > J c D ab } , (13) 

where the chemical potential E va \ < /i < m. Note that the order of D ab and J c is 
unambiguously dictated by the time ordering of the path-integral. One can always 
choose the basis \h > (by a transformation of \n >'s in each G sector separately) such, 
that: 

< n | I k >=< k \ T b j \n>, < h \ r c | k >= - <k \ r c \n> . (14) 



If at this point one does not pay attention to the order of D a b and J c in eq. fljjD , one 
can easily get a zero answer for (Aj)^ in view of the symmetry properties of eq . (|14D . 
However, keeping the right order of the collective operators and making use of the 
property [J a , Djb] = ie a b c Dj C , one finally arrives at the following expressions for the 
collective axial constants operators: 

si 3) = -(« + f) £33, gf=]h- (is) 

Note that the 1/N C correction corresponding to (3/1 has the same group-theoretical 



factor .D33 as the leading term in accordance with the general theorem of Ref. [19| The 
valence part of a entering eq.(|l5|) and following from eq.flTOD is given by: 

N c 

a va \ = — -(val I ar | val). (16) 
3 

Furthermore from eq.(fl3|): 

a N c ■ (val I (JiTj I n) (n I r fc I val) . 

/?vai = -— lUjk 2^ B B signE n , (17) 

n _, val h n - ^ va i 

N c ^ (val I a I n)(n \ r \ val) 

7val = -— 2^ 7T— p • (18) 

u n^val - C/n -^val 

Now we will proceed to define the constituent quark limit more precisely. For the 
sake of simplicity we will use a somewhat unphysical profile function P = 7i6(r — r) 
and then tune r — > 0. In this limit the energy of the valence level E wa \ — > m and the sea 
contributions to all quantities vanish. Simultaneously we will take the nonrelativistic 
limit, i.e. we will systematically neglect small components of the spinors (^) and (|5|). 
As we shall illustrate later the results are not changed for a realistic smooth profile 
function. 

For the 9 profile the term proportional to sinP in eq.(^) vanishes and the radial 
functions C(r) = D(r) = for allr, for boundary conditions i — 1. Similarly A(r) = 
B(r) = for i = 2. Therefore for given energy E n , G and M we have four independent 
solutions corresponding to two parities and two boundary conditions: 
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(x\n,G,M,+,l)=( '5 G "ffe!l' (x\n,G,M,+,2)= ( ^) 1> 

(* | n, G, M, - 1) = ( ) , (* | n, G, M, - 2) = ( 'j^Jfctf ) ■ (19) 

V AiG '— 'G,M / V ^nG '— 'G,M / 

For the valence level "val" corresponds to (n = 0,G = 0,M = 0,U = +,i = 1). Note 
that Sq~ ' ± ' ) identically vanish, so for G = there is no solution corresponding to the 
second boundary conditions i = 2. 

The nonrelativistic limit of eq.(||) can be obtained by approximating E w m. Then, 
for r > ro A + and C + become large components of natural parity spinors, whereas for 
unnatural parity large components are given by B~ and D~ . For r < r the upper and 
lower components are interchanged, since the sign of the mass term in eq. (P) flips, once 
r crosses ro- Therefore, the nonrelativistic limit corresponds to keeping B + = D + w 
for r > r and A~ = C~ ~ for r < r . The normalized solutions can be obtained if 
the system is confined in a box of radius R. Then the normalization condition for the 
valence level in the nonrelativistic limit reads: 



dr r 2 | B+ | 2 + dr r 2 | A+ Q \ 2 = 1. (20) 



r R 
J J.. .2 I n+ / .2 

r-o 

Obviously, in the limit ro — > the first integral has to vanish and the normalization 
is given entirely by the second term in eq.(^). One can similarly normalize other 
solutions corresponding to the nonvalence orbits. 

Now, in order to calculate constants a va i, /5 va i, 7 va i and the valence part of the 
moment of inertia we have to know the matrix elements of intrinsic operators between 
the states of eq.([L!|). Constant a va \ can be easily evaluated by means of the Wigner- 
Eckart theorem: 

«vai = ~f (- jdr r 2 | B+ | 2 +3 j dr r 2 | | 2 ] . (21) 

V r / 

In the limit r — > the first integral vanishes and a v3u \ = —N c , due to the normalization 
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condition (p0|). The proton-spin-up matrix element of the collective operator D33 is 
equal to —1/3, so the leading term of is equal to N c /3. 

In order to evaluate the next term in eq. fll5|) let us first observe that the matrix ele- 
ments (n |r|val) do not vanish only for unnatural parity G=l state n = (n, 1, M', — , 2) 
- this can be again easily verified by the Wigner-Eckart theorem. Here M' = 0, ±1 
depending on whether we take r 3 or r± respectively. Since the ratio /5 va i// va i is domi- 
nated by the lowest energy level which has the energy Eq ~ m the sums in ([18]) and 
(H) reduce actually only to the one term. After calculating ze^-fc(val| 0*7^1 n), where 
again n — (0,1, M r , — , 2), one gets finally: 

N 1 (1 R t 
(3 val = -f lim — / dr r 2 B+C 01 - / dr r 2 A+D c 

3 E ^m _fc n — t, j \ J J 



\0 r 

to R 



X 



J dr r 2 B+C 01 + 3jdr r 2 A+ D 01 . (22) 



\0 ro 



Similarly: 

ro R 



= -f \im m — _^ I rfr r^+Col - | rfr r 2 A+L>oi . (23) 



v0 r 

The constituent quark limit (CQL) consists in taking r — * and _E va i — > m: 

Pval CQL 



Ival 

Proceeding in a similar way one gets: 

N c . ,. 1 



-2. (24) 



ro -R 



lim — / dr r 2 B+ C 01 - f dr r 2 A+ D Q 

\0 r 
/ r R 

xljdr r 2 B+ Q C m + 3 J dr r 2 A+ D^ ) , (25) 



\0 r 

so that: 



Tval CQL 2^ ^ 
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Altogether we find: 



£)CQLiVc + 2 ^ ^( 0) cQL i ^ 



9a 3 



in perfect agreement with the CQM results. One can convince oneself by a direct 
numerical calculation, that this result holds for any realistic soliton profile, whose size 

(3) 

r — > 0. This is illustrated in Fig.l, where the g A dependence on r is plotted for 
P = 2 arctan((ro/r) 2 ). 

Now, in order to investigate the large soliton limit, it is enough to observe, that 
in this limit a ~ (3 [15[| in the sense of radial dependence and both scale like Tq [|15f 



whereas / scales like Tq ||20|| . Therefore for tq — > oo the new 0(1) terms are suppressed 



and the results of the Skyrme model are parametrically reproduced. However for large 
but finite tq, the correction exists and differs from the Skyrme model in contrast to 
the chiral bag calculations in Ref. [TJJ. In the case of g A the large soliton limit of 7 
goes as inverse powers of r [2~I|, so that it is vanishing and the Skyrme limit g A = is 



reproduced. So we have in fact a smooth interpolation between the CQM and Skyrme 
model. 

Summarizing: the axial coupling constants g^ and g A evaluated in the Nambu- 
Jona-Lasinio model in a small soliton limit agree with the Constituent Quark Model. 
Since the N JL model is a nontrivial dynamical model with valence quarks and polarized 
Dirac sea, and the CQM is based on purely kinematical assumptions concerning the 
symmetries of the wave functions, this agreement is by no means obvious. It is only 
possible because of the subtle cancellation of the overlap integrals of a priori very 
different operators. It is also important that the overlaps inside the soliton vanish 
with the soliton size approaching 0. Moreover, if the 0(1) corrections coming from the 
time- ordering of the collective operators were not taken into account, the agreement for 

(3) 

g A would never be possible. For it is hard to believe that other corrections, like meson 
loops for example, which require extra regularization and envisage completely different 
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dynamics, might bring about the similar cancellation (although this statement requires 
obviously quantitative analysis). Therefore we consider the results of this paper as a 
further confirmation of the correctness of the procedure presented in Refs. []13| — pT5[1 , 
where the detailed derivation of the 0(1) rotational corrections has been presented. 
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FIGURES 
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FIG. 1. The axial vector coupling constant g A in dependence of the size ro of the chiral 
field for the lowest order Q° and next to leading order 17° + S7 1 correction. The constituent 
quark mass is m = 370 MeV. 
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